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Abstract 



^ , The theory of self-consistent effective interactions in nuclei is extended for a 



system with a velocity dependent mean potential. By means of the field cou- 
pling method, we present a general prescription to derive effective interactions 
which are consistent with the mean potential. For a deformed system with 
the conventional pairing field, the velocity dependent effective interactions are 
derived as the multipole pairing interactions in doubly-stretched coordinates. 
They are applied to the microscopic analysis of the giant dipole resonances 
(GDR's) of ^^^'"^^^Sm, the first excited 2+ states of Sn isotopes and the first 
excited 3~ states of Mo isotopes. It is clarified that the interactions play cru- 
cial roles in describing the splitting and structure of GDR peaks, in restoring 



1 



the energy weighted sum rule, and in reducing the values of B{EX). 
PACS number(s): 21.30.Fe, 21.60.Ev, 24.30.Cz, 27.60.+j 



Typeset using REVT^ 



I. INTRODUCTION 



The understanding of nuclear collective excitations has been one of the most important 
subjects in the nuclear many body problem In microscopic analyses of such excitations, 
separable multipole interactions have been introduced and applied extensively for spherical, 
deformed and rotating nuclei. The separable multipole-multipole interactions originate from 
long-range correlations in particle-hole channels, and the physical meaning of these interac- 
tions has been clarified by Mottelson in terms of the core polarization phenomena 0. On 
the other hand, the multipole pairing interactions originate from short-range correlations 
in particle-particle channels, and there are many works concerning the physical meaning of 
these interactions These interactions have been widely used and have been playing 

crucial roles in the study of nuclear structure, but theoretical foundations for the origin of 
such effective interactions in both the particle-hole and particle-particle channels have not 
been established enough from a unified physical picture. In particular, we consider that it 
is still an open and interesting problem to determine the proper form of the deformation 
dependence of the pairing interaction. It is the main purpose of the present paper to report 
a general prescription to derive effective interactions which are consistent with the mean 
potential, and we will present a unified derivation of the multipole-multipole interactions in 
particle-hole channels and the multipole pairing interactions in particle-particle channels on 
an equal footing (preliminary reports on this subject can be seen in Refs. [||,||). 

Since a nucleus can be regarded as a spatially and energetically saturated self-sustained 
system with a relatively sharp boundary, one can assume that the following condition, which 
will be called as nuclear self- consistency is satisfied quite accurately: The shape of 

the mean potential and that of the density are the same even when the system undergoes 
collective motions. The concept of nuclear self-consistency, which is much more stringent 
than the Hartree self- consistency, has played important roles not only in describing an 



equilibrium nuclear shape p2| , pr3| and in deriving effective interactions for a spherical system 
|1|J3], but also in deriving effective interactions for a deformed system |TO|Tl|], and higher 
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order effective interactions [^,0,0. These interactions have been apphed successfully 



to the calculations of the properties of the low-lying vibrational states and high-frequency 
giant resonances, etc. [ p!0| , pJ] ,p!6|-p3| . It is found that the effective interactions derived by 
the rigorous application of the nuclear self-consistency are much more reliable than the 
conventional multipole-multipole interactions |I|,0,|I^,|TT|] . 

However, it has been pointed out that the spurious velocity dependence in a single 
particle mean potential, which violates the Galilean invariance of the system, should be 
removed It is well recognized that the velocity dependent terms affect the mass 

parameters of collective motions and thereby such quantities as the distribution of transition 
strengths and the absolute value of the energy- weighted sum-rule, etc., become unreliable 
Plj3|j5H7|, |2^p5| . The realistic nuclear potential in fact contains velocity dependent terms 



such as the BCS pairing field, / ■ s and P terms in the Nilsson model, etc. Some parts 
of them are responsible for violating the Galilean invariance and the classical sum-rule, 
and therefore these symmetries have to be restored properly. From this point of view, 
Pyatov et al. developed a simple and powerful method to introduce additional interactions 
which restore such broken symmetries in RPA order 0,^. Since then, symmetry restoring 
effective interactions have been studied by several authors |[5|-|7|P^|-pU[ . In this paper, by 
investigating the coupling between collective and single-particle degrees of freedom, and by 
using the concept of nuclear self- consistency and local Galilean invariance of the system as 
important guiding principles, we will present a systematic method to derive self-consistent 
effective interactions in nuclei. 

In sect.||, the field coupling method ||l|,pT| is applied to a spherical system with a velocity 



dependent potential. In this framework, we derive self-consistent velocity dependent effective 
interactions by estimating the coupling between the collective displacement of nucleons and 
the mean field. In sect.|Tl|, we extend the field coupling method to a deformed system 
and present a simple derivation of the doubly-stretched multipole-multipole and multipole- 
pairing interactions. In sect.fV], we investigate some fundamental properties of these self- 
consistent effective interactions. It is shown to be essential to express effective interactions in 



doubly-stretched coordinates for restoring some broken symmetries and also for the natural 
description of GDR's in deformed nuclei. In sect.jV], we report the results of numerical 
calculations in RPA of the GDR's of ^^^'^^^Sm, the first excited 2+ states of Sn isotopes and 
the first excited 3~ states of Mo isotopes, by using the self-consistent velocity dependent 
effective interactions. 



II. GENERAL FRAMEWORK OF FIELD COUPLING METHOD 

A. Field coupling between collective distortion and mean potential 

The method to study collective motion that arises from the action of the field coupling 
in a system with a degenerate one-particle excitation has been well developed by Bohr 
and Mottelson It provides a self-consistent method to construct the relevant effective 
interaction by identifying the field coupling as the Hartree field of the interaction. Based 
on the method, which will be referred to as the field coupling method, the theory of self- 



consistent effective interactions in nuclei has been developed [ 10, 11, 14 1. As a result, it has 
been shown for example that the conventional multipole interaction model must be improved 
to satisfy the nuclear self-consistency in deformed nuclei, resulting in the doubly-stretched 
multipole interaction model. Here we will extensively apply the method to a system with the 
velocity dependent mean potential. Some examples of the velocity dependent field coupling 
were discussed in Ref. [|] in connection with the analysis of the high-frequency quadrupole 
modes and of the center-of-mass mode. 

In a nuclear system, particle motions and collective motions are essentially coupled with 
each other in order to achieve self-consistency between the potential and the density distri- 
bution. If a collective mode excites in a nucleus, the positions and momenta of nucleons 
in the system are displaced by the mode accordingly. Then the corresponding change in 
the nuclear density distribution gives rise to a violation of the self-consistency settled before 
the displacements are switched on. In this way additional field couplings are induced in the 
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system in order to restore the nuclear self-consistency. 

Let us start by considering a 2'^-pole collective shape oscillation mode of a spherical 
nucleus. We assume that this mode is characterized by the collective displacement of the 
nucleonic field variable as 

r ^ r + 5r, 5r = - J] a^^Vg^M (A > 1), (1) 

where are the collective amplitudes. The collective velocity field associated with this 
mode is expressed as 

v{r) = -Y.al^VQ^,, (2) 

which is irrotational and incompressible. We will require that nucleonic velocities enter- 
ing into the velocity-dependent single-particle potential are to be measured relative to the 
collective velocity v{r) so that the potential is to be invariant under the local Galilean 
transformation 

p ^ 'P + Sp, Sp = Mv{f). (3) 

Then the variation of the average one-body potential 5V and that of the density distribution 
5p produced by the oscillation are determined from the conditions 

V{f+ 6f,p + 5p) = Vo{f,p}, 

p{r + Sr,p + Sp) = po{r,p), (4) 

where Vo{r,p) and po{r,p) are the potential and the density at the original equilibrium point 
in the phase space and are assumed to be spherical, while V{f,p) and p{r,p) include the 
effect of the oscillation. In this paper we assume that these quantities expressed in the phase 
space are well defined by using some appropriate semiclassical method such as the Wigner 



transformation 31 



The conditions of Eq.(^ provide relations between 5V and 5p through the displacement 
vectors 5r and 5p as 
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V{f,p) = Vo{r,p) + SV{f,p), 

p(f,p) = po{r,p) + 6p{f,p), (5) 



with 



6V{f,p) = 6Vr + 6Vp = -6r- Wo - 6p- VpVo, 

5p{f,p) = 5pr + 5pp = -5f- Vpo - Sp- Vppo, (6) 



where notations are defined by 



\dx^ dy^ dz J ^ ^ \dpx^ dpy' dpz 



In Eq.(l^), for simphcity, only the leading order terms in and are retained assuming 
that the vibrational amplitudes are small. Inclusion of the non-linear field coupling coming 
from higher-order terms in has been performed in Ref. resulting in higher-order 
(many-body) effective interactions. The 5Vr and 5Vp can be expressed in the standard form 
of the field couphng [jl| as 

5Vr = KxYl «v^am' = — VQam ■ VV^o, (8) 

6Vp = k,J2 «I/am, K = ^^Qx, ■ Vp^o. (9) 
^l ^\ 

Now we consider that the field couplings 5Vr and 5Vp are produced as the Hartree field 

of two-body interactions 

H^nt = f E FlF,, + f E FlF,, ^ Ht' + iff, (10) 

and introduce self-consistency conditions as 

al^ = A< Fl >, al, = A< fI > . (11) 

Here the average of a one-body operator with respect to the modified ground state | > 
corresponding to the density p is calculated as 



A 

A<F >=< ^F(i) >= / Fp{f,p)d^rd^p, (12) 
i=i •' 

while that with respect to the original ground state |0 > corresponding to the density po is 
given by 

A 



A<F >o=< >o= / Fpoir, p)d\d^p. (13) 

i=i 

Equation (|TTp can be calculated as 

«Im = / ^L^Pr d'^rd^P 

= -al^A < (VFl . VQ,,) >o, (14) 
= Mal^ J FlVQx, ■ Vppo d\d'p 

= -Mal,A < {V,Fl . Vg^) >o, (15) 
where we have assumed the time reversal invariance of the density po which guarantees 

J F^Spp d\d'p cx A < V,Fl ■ VQam >o= 0, (16) 

J Fl6pr d'rd'p oc A < VFl ■ VQx, >o= 0. (17) 



We thus obtain 



Kx = -A<V (VQi^ ■ VVo) ■ VQam >o, (18) 



-M^A < Vp (VQl^ ■ V,Vo) ■ VQam >o • (19) 



Since Qx^ does not depend on the momentum p, Fx^ of Eq.(|^) and kx of Eq.(|T9D can be 
expressed in terms of the Poisson bracket as 

^Am = ^^{Qam,K)}, (20) 



Rx = -M' < {Qx„ {Ql^, Vo}} >o . (21) 
In the above treatment, we have assumed the potential Vq can be expressed as a sum 

A 

over individual particles: Vq = J2 ^oi'^^ii Pi)- However the monopole pairing potential 

i=l 

v,a.r = -^(4 + ^o), 4 = E 44 = 2 E 44 (22) 

is not expressed in this form. Here cl^ creates a nucleon in the state 7, 7 is a time-reversed 
state of 7, and A is a gap parameter. For such a potential we will replace the Poisson 
brackets in Eqs. (pOD and (^) by the corresponding commutation relations as 

M 



^[Qx,,Vo], (23) 



(-^ <[Qx„[Vo,Ql]]>o. (24) 



Then the corresponding part of the effective interaction becomes 

This type of separable effective interactions have been discussed by several authors in con- 
nection with the symmetry- restoring treatment of the nuclear Hamiltonian ||5|-|7,26- 0O|] . 



B. Two simple examples of self-consistent effective interactions 

It is worthwhile to briefly review the derivation of the self-consistent effective interactions 
in the case of a simple mean potential such as a harmonic oscillator potential |Tl| or a pairing 
potential [0. 

1. Spherical harmonic oscillator potential 



Let us consider a simple situation that the equilibrium potential Vq is given by the 
spherical harmonic oscillator potential 
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Vo = VUr) = \Muoy. (26) 

In this case, we do not have to consider the velocity dependent field coupling since the 
potential does not have any velocity dependence. From Eqs.(||) and ([18|) we obtain 

Fx^ = —Mul f ■ VQam = —XMul Qxf,, (27) 

Kx = -XMujIA < VQ{^ ■ VQx^ >o 

= _Ma;2^!i^A±llA < r2^-2 (28) 

4:71 

where we used the relations 

HQlQx.) = '-^r^'-\ (29) 

Thus the self-consistent effective interaction coincides with the conventional multipole- 
multipole interaction 



2' . 



with the self-consistent strength 



Air Muj, 



2 



self ^ ^y^^O (^-.s 



2. Pairing potential 
Let us introduce a one-body operator 

= J2< a\0\P > cicp, (32) 

a/3 

whose time-reversal property is assumed to be 

f of -1 = {-fO^ (33) 
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or equivalently 

< a|0|/3 >=< fa\0\f(3 >= {-f < l3\0\a >, (34) 

where T denotes the time-reversal operator, and (— is a short-hand notation for the time- 
reversal phase and is either +1 or —1. The commutation relation between the operator O 
and the pairing field Vpair is given by 

[0,Vpa^r] = A E < "l^l/^ > (c^^ -C^C^). (35) 

^ a/3 

From this relation, it is easily recognized that the pairing field satisfies the translational 
invariance but violates the Galilean invariance: 

[P, VJ,a^r] = 0, (36) 

[R, Vpa^r] = A ^ < c^W > (c^^ " ^/gCa) (37) 

with 

^ = ER;=EH^V,), i?=^Er-;-. (38) 

More generally, if O is an operator depending on coordinate variables, then Eg . (p5D means 
that the monopole pairing field is not invariant under a local Galilean transformation. In 
this sense the pairing field is considered to be a velocity dependent mean potential [0] . For 
a system with such a velocity dependent potential, there arises a velocity dependent field 
coupling as is discussed in subsect.fT^. 

Let us investigate the interaction of Eq.(^3|) in more detail. By substituting Vpair into 
Vo and by identifying the original ground state |0 > with the BCS vacuum state, the basic 
quantities in Eqs. (P^D and (^) are expressed as 

[Qx,,V,an] = -A{Pl- P~\ (39) 



< [Qxfj. , [Vpair, Qx^] ] >0 



A^E(^ + ^)l<"l^?A.|/3>r, (40) 
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where 

4. = E<«IQA,|/5>ctct, (41) 

and we have used the famihar relation between the coefficients of Bogohubov-Valatin trans- 
formation and the quasiparticle energy: UaVa = A/2i?Q. Thus the velocity dependent 
effective interaction of Eq. (P5|) can be obtained as 



= -lof' UPl - PtJ(P>^. - Pi) (42) 

with 

Gf' /T.\i^ + ^)\< «l^?A,|/3 > r • (43) 

This interaction is a natural extension of the dipole-pairing interaction obtained by Bohr 
and Mottelson |l[] and Pyatov and Salamov [^] to a general 2'^-pole mode. Equation (^Sf ) 



has the same structure as that of the Belyaev identity Further discussions for this 

interaction can be found in Refs. ||7|,|20[ . 



III. SELF-CONSISTENT EFFECTIVE INTERACTIONS IN DEFORMED 

NUCLEI 

Now, let us apply the procedure of sect.|| to a deformed system. For simplicity, we 
assume that the main part of the equilibrium potential is described as a deformed harmonic 
oscillator with frequencies Ux, ojy and Uz. In terms of the doubly-stretched coordinates 
IO| , p!l|] defined by 



X = — X, y = —y, z = — z, (44) 

UJq Uq LOq 



the deformed harmonic oscillator potential can be expressed in a spherical form as 



2 2 

UJ,,Z 

^Mu'^irr = VUr"). (45) 
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As one of the possible and plausible collective shape oscillation mode in the deformed 
system, we will introduce a doubly-stretched 2'^-pole mode characterized by 

r"^r" + 5T\ 5t" = -Y.^I^^"Q"x, (A > 1), (46) 

with 

rather than the conventional 2'^-pole mode characterized by Eq.(|l]). In fact, as discussed 



in Ref. [TTI based on the Thomas- Fermi theory, the doubly-stretched mode represents great 
improvements over the conventional mode in the sense that it satisfies (i) the constancy of 
the Fermi energy (which is equivalent to the saturation condition), (ii) the separation of the 
center-of-mass motion, (iii) the condition for a fluctuation around the deformed equilibrium 
shape, and (iv) the self-consistency between the nucleonic density and the potential, etc. 
The collective velocity field associated with this mode is expressed as 

v"{n = -i:^i,^"Qi,^ (48) 

which is a natural extension of the irrotational and incompressible flow of Eq.(0). 

Now we will extensively apply the requirement of the local Galilean invariance to the 
deformed system. Namely, we will require that in terms of the doubly-stretched coordinates 
nucleonic velocities entering into the velocity-dependent single-particle potential are to be 
measured relative to the collective velocity v"{f") so that the potential is to be invariant 
under the transformation 

p" ^ p" + Sp", 5p" = Mv"{f"), (49) 

where the doubly-stretched momentums are introduced as 

Px = — Px, Py = —Py, Pz = —Pz- (50) 

Ux UJy UJz 

Then by use of the relations d^r = d?r" and d^p = d^p", which are equivalent to the volume 
conservation condition uOxUjyUJz = uJq, we can follow the similar procedure as given in sect. HI 
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to derive self-consistent effective interactions in deformed nuclei. To do this, in the present 
case, we should understand that all the r, p, V and Vp appearing in sect.^ are to be replaced 
by the corresponding r ", p ", V" and with 

"^'^^{wi'WfWi)' ^^^^ 

For example, the field coupling induced by the collective displacements of Eqs. (|46|) and (^91) 
becomes 

Vif",p") = Voif",p") + 5V{f",p") = V^' + 5V" (52) 

with 

5v" = 5v;' + sv;; 

= E E (53) 

K = -^"Ql ■ V 'X, (54) 
K = -^"Ql ■ V = -{Ql, K}- (55) 

From the self-consistency conditions 

ax^ = A< i^; >, ax^ = A< F'^^ >, (56) 
the coupling strengths are determined as 

^x = -A< V" {V"QIJ ■ V'X) ■ ^"Ql >o, (57) 

= -M'<{Ql,{QVX}}>o, (58) 
and the self-consistent effective interaction to be used in the deformed nucleus is given as 

TT _ '^-^ t 17" I '^^ \^ TP" t TP" 

^Ht^ + H^\ (59) 
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For the case of Vq = Vho{r"), the effective interaction becomes 

1 

with 



H^x^ = ~xf'j:QVQl (60) 



^ff = -± (61) 

2A + 1A< (r")2^-2 >o' 

and in the presence of the pairing field an additional effective interaction is derived as 



H? = -Icf' E(nV - - P^J) (62) 

with 

= 1 /E + ^)l < c^\Ql\/3 > r • (63) 

These interactions play crucial roles in restoring broken symmetries of the system in RPA 
order, the detail of which will be discussed in the next section. 



IV. FUNDAMENTAL PROPERTIES OF SELF-CONSISTENT EFFECTIVE 

INTERACTIONS 

A. Translational invariance of a deformed system 

The consistency of the residual interaction with the shell model potential and the method 
to restore the translational invariance of a nuclear many-body system was discussed by 
Pyatov and Salamov |^ for the case of the spherical oscillator potential. Here we will briefly 
examine this problem in a deformed system. 

Let us consider the oscillator Hamiltonian for the deformed nucleus 

n«-t{i} + yur';)}. (64) 

Since Hq is a local one-body Hamiltonian, it breaks the translational invariance of the 
system. In fact, we have 
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[Ho, P"] = ihAMulR" (65) 

with 

p" = E p/' = E (-^ w;) , ^" = ^ E ^/ ■ (66) 

i=i i=i i=i 

To recover the translational invariance, we need a counter term which cancels the above 

commutation relation. For a spherical nucleus, as is well known, such a counter term comes 

from the conventional dipole interaction. Now we will show that for the deformed nucleus, 

the doubly-stretched dipole interaction plays the role to recover the translational invariance 

of the system. 

First of all, for the doubly-stretched dipole interaction 

^A=i = -fE(^i«-Qi'(i)), (67) 



we can verify 



[V^=u P"] = -ihAMu'o^R". (68) 

Xl 

Then the total Hamiltonian H — Hq + Vx=i satisfies 

[H, P"] = ihAMul (l - j R". (69) 

Therefore if the strength Xi of the doubly-stretched dipole interaction is set equal to its 
self-consistent value xT^^ — 4:'tMu!q/3A, we obtain [H,P"] — exactly. In this case, the 
total Hamiltonian can be expressed as 

1=1 

which explicitly guarantees the translational invariance of the deformed system. 



B. Restoration of the local Galilean invariance 

As is discussed in the previous sections, generally a phenomenological potential, com- 
prising velocity dependent terms such as the 1^ term, the I ■ s terms and the pairing field, 
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etc., does not commute with arbitrary coordinate operators. If we chose the operator to be 
the multipole operator, such a situation is expressed from Eq.p3|) as 

M 



[Qv'W,P)] = ^^;^o, (71) 



which means that the potential Vq violates the local Galilean invariance under the collective 
multipole oscillation. Here the doubly-stretched operators are used for deformed nuclei. For 
spherical nuclei, of course we can omit the double primes. 

Now along the line of the general method of restoring the broken symmetry ||5|j2^ , we will 



show that the self-consistent velocity dependent effective interaction, i.e., of Eq.(^), 
plays the role to restore the local Galilean invariance of the system under the random phase 
approximation (RPA). In the RPA order we can verify 

, H^^\pA = I {[Qi , f'^jupa, ^;}^ , (72) 

with 

[Ql , fVUpa = ^ < [Ql , [QV > ^o]] >o, = (73) 



Thus we obtain 



[Q'i,,Vo + Hi'^UpA = 0. (74) 



C. Simple model analysis of GDR of normal nuclei 

It is worthwhile to point out that the doubly-stretched interaction model is powerful 
and plausible also for the description of some iso- vector modes. In confirmation of it, let 
us briefly review the simple model analysis of the splitting of GDR in an axially symmetric 
deformed nucleus The model Hamiltonian is assumed to be 

H = Ho + vZT'\ (75) 

where Hq is the deformed oscillator Hamiltonian of Eq.(|5^ while is a residual iso- 

vector dipole interaction. Here we will parametrize the shape of the nuclear potential as 
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oj^ = = Uy = ujQ{e)(l + e/3), = ujo{e){l - 2e/3) (76) 

with 

For comparison we will introduce two types of iso- vector dipole interaction, one is an ordinary 
type and the other is a doubly-stretched type. The interaction in the latter case is given by 

yZ:'^ = E x^^K'\Q';Kr.y {q'Ikt.) (78) 

^ K 

and we will parametrize the force strength as Xix"^^ — "^xT^^ i where xT'^^ is the self- 
consistent strength of the iso-scalar dipole interaction. The typical value of ^ estimated 
from the symmetry energy term in the mass formula under the Fermi gas approximation is 
about 3 111]. 

Under the RPA, the excitation energy of each K-component of the GDR is obtained 
analytically as 

^11 = ^GDR I 1 + 



3(1+0. 
2e 

3(1 + 0, 



^10 = ^GDR (1 - 777^^ I ' (79) 



for the ordinary interaction, while 

^^11 = ^GDR (1 + e/3) , = nGDR (1 - 2e/3) , (80) 

for the doubly-stretched interaction. Here ^Igdr is the resonance energy of a spherical 
nucleus given by 



^GDR = V 1 + ^^0, (81) 

which, in both cases, is compatible with the experimental systematics of ^gdr ~ 80 A~^^^ 
[MeV] if we put ^ = 3. For the doubly-stretched interaction, the total energy splitting 
between the K=0 and K=l components of the GDR is e in units of ^Igdr and is independent 



of ^, which is consistent with the simple classical geometrical relation of {uj± — uJz)/ujo = ^• 



This is in good agreement with the systematics of the experimental observation p5|-|37[| . On 
the other hand, for the ordinary interaction, the splitting is too small by a factor of 4 for 
^ = 3. Thus the doubly-stretched interaction model seems much more improved than the 
ordinary one also for the iso-vector dipole mode. 

D. Simple model analysis of GDR of superconductive nuclei 

Effects of the inclusion of the pairing correlation on the properties of giant resonances 
of superconductive nuclei have been studied by several authors |||,|^,|^, and the shifts 
in excitation energies, the changes in energy- weighted sum rule (EWSR), etc., have been 



observed. However as indicated in Refs. [24,25], some of them seem to be spurious due to 



the violation of the Galilean invariance of the system. Here, by use of a simple schematic 
model, we will verify the effects of the pairing correlations on the structure of GDR, and will 
show that such spurious effects can be remedied by including the dipole pairing interaction. 
The model Hamiltonian is assumed to be 

H = Ho + vZT''' + E {^Pa^r ' + H^^l,} ^ , (82) 

T 

where Hq is the oscillator Hamiltonian of Eq.(|^), v!^'^'^^^ is the iso-vector dipole interaction 
of Eq.(^), Vpair is the monopole pairing potential of Eq.(|^), and H^l^ is the dipole pairing 
interaction of Eq.(|62D. Here and in the following, the summation index r is to be taken 
over the proton and the neutron. The force strength of the iso-vector dipole interaction is 
fixed as ^ = 3, while that of the dipole pairing interaction is set to be its self-consistent 
value when it is included. We perform quasiparticle RPA calculations for a schematic model 
system with N=Z=20. It must be noticed that the purpose here is not to compare with 
the experimental data of ^°Ca but to investigate the fundamental properties of the dipole 
pairing interaction from the purely theoretical point of view. 

Table | shows our results on the electric dipole strength distributions calculated by 
assuming the system to be spherical. To study the effect of the pairing correlation, we 
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artificially change the value of A. For the simplest case of A = (a), the GDR is located at 
24.0 MeV, and the energy-weighted El transition strength of this state exhausts its classical 
sum rule value of 

S{El),i,,, = — — ^e\ (83) 
V >class 4^2M A ^ ' 

When the pairing gap is set to be A = 1.0 MeV (b), reflecting the situation that the 
quasiparticle states are no more degenerate, the GDR splits in spite of the assumption that 
the system is spherical. There are mainly three components of the GDR ( 21.4, 23.6 and 
27.5 MeV). The center of the GDR shifts upward about 1 MeV compared to the case of 
A = 0, and the EWSR is overestimated by about 10% relative to the classical one. In order 
to remedy this situation, we switch on the dipole pairing interaction (c). In this case, the El 
transition strength concentrates again on a single state at about 24.0 MeV, and the EWSR 
recovers to its classical value. The fact that the pairing correlation destroys the order in the 
structure of the GDR can further be emphasized by increasing A to 2.0 MeV, though the 
value is not realistic. In this hypothetical situation (d), the transition strength splits into 
mainly four components at 20.8, 23.6, 30.0 and 51.3 MeV, and the EWSR is overestimated 
by about 36%. Even in such an extreme situation, if we additionally include the dipole 
pairing interaction (e), the spurious effect of the pairing interaction can be removed and the 
El strength is concentrated on a single state at 23.9 MeV to recover the EWSR. 

Figure |l] shows the corresponding results obtained by assuming that the model system 
is deformed to be axially symmetric shape oi (3 = 0.4. The continuous strength function, 
representing the transition strength per unit energy, is constructed by using the Lorentzian 
weight function 

P{U)-UJ^) = -— I r2 2 - (^^) 

TT [UJ^ — Ujf^Y + i 

In the present model calculations, we choose the width to be F = 1.0 MeV only for the 
sake of not wiping out the fine structure of the resonance. Because of the deformation 
of the system, the GDR splits into K=0 and K=l components. For the case of (a), the 
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GDR shows rather comphcated structure. The reason for it can be traced to the violation 
of the Gahlean invariance for the pairing field and we can eliminate such a spurious effect 
by restoring the broken invariance. In fact, if we additionally include the dipole pairing 
interaction, the structure of the resonance becomes simpler both for K=0 and K—1 modes 
(b). Furthermore, the total energy sphtting between the K=0 and K—1 modes, relative to 
the average energy of these modes, becomes approximately equal to /? which is consistent 
to the simple classical geometrical relation explained before. 

V. NUMERICAL RESULTS 

In this section, we report some characteristic results of numerical calculations in RPA 
of the GDR's of ^^^'^^^Sm, the first excited 2+ states of Sn isotopes and the first excited 3~ 
states of Mo isotopes. It should be noted here that our present calculations contain essen- 
tially only one free parameter in the following sense; the single particle bases are constructed 
from the Nilsson + BCS model with standard parameters; the strengths of the velocity de- 
pendent effective interactions such as the multipole pairing interactions are fixed to be their 
self-consistent values when they are included; only the strengths of the multipole-multipole 
interactions are adjusted under the condition that a common value of ^, X2 and Xs is to be 
adopted for all the isotopes of Sm, Sn and Mo, respectively. 

Since the positions of band-head states are very sensitive to the choice of the single 
particle energies as are generally observed in the RPA calculation of vibrational states in 
deformed nuclei, one can improve the fit to experiments by adjusting the single-particle 
energies, and further improvement can be obtained by slightly varying the strengths of 
the effective interactions around the vicinity of the predicted self-consistent values for each 
isotope. The fit to experiment we obtained is insufficient and an improved fit could have 
been obtained if these parameters were treated as adjustable as well. We have not done 
so in the present work, because the purpose of the present numerical investigation is to see 
first to what extent our theory works without much playing around with parameters and to 
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provide understandings rather than precision tools for the fitting of experimental data. 



A. El strength distributions of J^4s,i54gj^ 

We here present the results of realistic calculations on the El strength distributions of 
i48,i54gj^ in the quasiparticle RPA. The model Hamiltonian is assumed to be 

H = Ho + Vmi + vZT'^ + J2 {^P-r -XN + ffg,}^ , (85) 

T 

A 

Vmi = E [-^^ ^0 {2(r- s) + /i(P- < ^"^ >)}]^ , (86) 
1=1 * 

which is essentially same as Eq.(^) except that the Nilsson potential Vnu is additionally 

included. 

We will study the effects of two kinds of velocity-dependent interactions to restore the 
Galilean invariance of the system; the one arising from Vpair and the other from the velocity 
dependent part of Vnu. The former, the dipole pairing interaction H^}^, is given by Eq.(|6^), 
while the latter, Hl^f \ is given by 

M=f = ^E'^i^U. (87) 

with 

Vl,= [Ql,VN^l], (88) 

K^ = -l/ <[Ql,[Vmi,Q';J]]>o. (89) 

Although it is known that there exists a term associated with the coordinate distortion of 
the spin-orbit potential given by Eq.(6-70) of Ref. we will not consider its effects in order 
to concentrate our present analysis on the effects of the self-consistent velocity dependent 
interactions. We will also neglect iso-vector corrections coming from the Nilsson potential 
to the iso-vector dipole-dipole interaction because of the same reason. 

22 



To fix the parameter ^ for the strength of the iso-vector dipole-dipole interaction in this 
mass region, we first calculate the case of ^^^Sm by assuming its shape to be spherical. The 
resultant value of = 2.7 is adopted also for ^^^Sm whose quadrupole deformation parameter 
is assumed and fixed as /5 = 0.35 in the following calculations. All other strengths of the 
interactions are fixed just as the self-consistent values which restore the Galilean invariance 
of the system. It should be noticed here that the doubly-stretched interactions are used 
for the deformed nucleus of ^^"^Sm. For the model space, we retain all the Nilsson single 
particle states with 2 < Nose < 7 for protons and 3 < Nose < 8 for neutrons. The Nilsson 



parameters are taken from Ref. ||13|- By using experimental binding energies of Ref. [38 



the pairing gap parameters are determined from the even-odd mass differences as A„ = 1.01 
MeV, Ap = 1.36 MeV for ^^^Sm and A„ = 1.07 MeV, Ap = 0.86 MeV for ^^^Sm. We use 
the Lorentzian distribution of Eq.(|84D to reproduce the resonance width. We choose F (in 
units of MeV) as 5.10 for ^^^Sm, and 3.25, 5.25 for the K=0, 1 modes of ^^'^Sm, respectively. 

Figures ^ and 0b show the El strength functions for ^^^Sm and ^^^Sm, respectively, 
taking the model Hamiltonian of Eq.(|55D but switching off the dipole-pairing interaction. 
Here we see that even if we use the doubly-stretched interaction, the calculated splitting 
between K=0 and K=l resonances of ^^^Sm is too small. Furthermore, as can be seen from 
the column a of Table the EWSR values exceed the classical values for both nuclei. As 
stated repeatedly, these difficulties stem from the mixture of spurious states arising from the 
broken Galilean invariance of Vpair and Vnu. In the following we show the results obtained 
by restoring the broken symmetry in two step. 

First, we study the effect of H^l^. The column b of Table || shows that the sum rule 
values approach the classical values. However, from Fig.|^, we see that the centers of the 
resonances shift to lower excitation energies and the structure in the lower peak region reveal 
unnatural shape. Finally we take into account H^l^ so that the Galilean invariance of the 
Hamiltonian is restored. As can be seen from Fig.^ and column c of Table 0, the unnatural 
resonance structure disappears and the EWSR values keep close to the classical limit. We 
note that the final results of Figs.^ and agree quite well with the Lorentzian distribution 
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of Figs.^ and ^ which fit the experimental data of photoneutron cross section 
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B. The first excited 2+ states in Sn isotopes 

We now consider the effects of the self-consistent quadrupole pairing interaction on the 
excitation energies -^(2^'') and the E2 transition probabilities B(E2) of Sn isotopes within 
the quasiparticle RPA. The model Hamiltonian is assumed to be 

H = Ho + Vmi + 14=2 + E {^P- - AiV + H^^},}^ , (90) 

r 

where Va=2 and -^1=2 quadrupole-quadrupole interaction and the quadrupole pairing 

interaction, respectively. 

The single particle model space is spanned by all the Nilsson states with 2 < Nose < 6 



for protons and 2 < Nose < 7 for neutrons. The Nilsson parameters are taken from Ref. [jT3 
and the deformation is set to be zero. To investigate the effect of the quadrupole pairing 
force on the B{E2; Og 2i) value, the strength of the quadrupole pairing force, G2, is fixed 
to be its self-consistent value when it is included, while that of the quadrupole interaction, 
X2, is used as an adjustable parameter to reproduce experimental E{2f) of 116,118,120^^ 



First, we adopt the pairing gaps determined from the even-odd mass differences. The 
adopted value of X2 for the best fit obtained without (with) the inclusion of the quadrupole- 
pairing interaction, H^x=2i 0.92 (0.90) in units of xr^'^- The calculated E{2i) and B(E2) 
are plotted in Figs.^ and P?, respectively, and corresponding numerical values are given 
in column a of Table |T|. These results show that H^x=2 necessary and important in 
reproducing experimental E{2i) and B(E2) values simultaneously. 

Second, we have also performed numerical calculations by fixing the proton energy gap 
to be zero. This is because the proton shell is closed in Sn. In this case, the adopted value of 
X2 for the best fit obtained without (with) the inclusion of iy{=2 0.940 (0.965) in units of 
■ The calculated values of E{2i) and B(E2) are plotted in Figs.||a and respectively, 
and corresponding numerical values are given in column b of Table |ITT|, where we see better 
agreements with experimental data compared to the case with Ap 7^ 0. 
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C. The first Excited 3 States in Mo Isotopes 

Here we will study the effects of the self-consistent octupole pairing interaction, H^l^, 
in the quasiparticle RPA calculations of the excitation energies E{3i) and the E3 transition 
probabilities of the first excited 3~ states in Mo isotopes. The model Hamiltonian is assumed 
to be 

H = Ho + V^u + Vx=3 + E {^P- - AiV + H^^l,}^ , (91) 

r 

where V\=s and H^ls are the octupole-octupole interaction and the octupole pairing inter- 
action, respectively. 

The single particle model space is spanned by all the Nilsson states with < Nose < 7 for 
protons and < Nose < 9 for neutrons. Energy gaps are determined from the experimental 
even-odd mass differences. To investigate the effect of the octupole pairing force on the 
B{E3;0g 3i) value, the strength of the octupole pairing force, G3, is fixed to be its 
self-consistent value when it is included, while that of the octupole interaction, X3, is used 
as an adjustable parameter to reproduce the experimental data of E{3i) in 94,96,98, 100 j^^^ 
The adopted value of Xs for the best fit obtained without (with) the inclusion of the 



octupole-pairing interaction, -^{=3, is 0.95 (0.97) in units of xT'^^ ■ 

The calculated E{3i) and B(E3) from the ground state to the 3j~ state are plotted in 
Figs.^a and |^, respectively, and corresponding numerical values are given in Table P^ . 
Here we see that the calculated B(E3) values are reduced and improved systematically by 
the effects of H^x=zi although the fit to experimental data is insufficient. 

VI. SUMMARY 

We have extensively applied the prescription to derive self-consistent effective interac- 
tions, needed for the unified description of the collective motions of atomic nuclei, especially 
to the system with velocity dependent mean potentials. As the guiding principles for this 
purpose we have imposed the conditions of nuclear self-consistency |llO| , |TT| and local Galilean 
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invariance |]l],^,|5|-0 of the system, implemented with the simple and transparent field cou- 
pling method developed by Bohr and Mottelson [^. 

The nuclear self-consistency requires that the shape of the mean potential and that of 
the density are the same even when the system undergoes collective motions, while the 
local Galilean invariance requires that the nucleonic velocities entering into the velocity- 
dependent single-particle potential should be expressed relative to the local collective fiow. 
In the field coupling method, the coupling between the particle motion and the collective 
field is identified, within the Hartree approximation, as the averaged one-body field of the 
effective interaction which we look for, and the coupling strength of it can be fixed by the 
above conditions. 

For the multipole collective shape oscillation modes in the harmonic oscillator potential 
with the monopole pairing correlation, we have derived the multipole-multipole interactions 
for particle-hole channel and the multipole-pairing interactions for particle-particle channel 
from the unified physical picture. In the case of deformed nuclei, it is shown that these in- 
teractions must be expressed in terms of the doubly-stretched coordinates so as to guarantee 
the conditions of nuclear self-consistency and local Galilean invariance of the system. 

The origin of the doubly-stretched multipole-multipole interactions have already been 
clarified and they have found many successful applications [|l0 i |Tl i [l6H23[| , while the origin 



of the doubly-stretched multipole pairing interactions are clarified in this paper on the 
same footing. We have applied the doubly-stretched multipole pairing interactions to the 
analyses of some collective states in Sm, Sn and Mo isotopes by means of RPA, and for 
the dipole mode we have also tested the velocity dependent effective interaction arising 
from the Nilsson potential. We have seen the effects of such velocity dependent effective 
interactions in the recovery of the classical El sum rule for GDR's of ^'^^'^^'^S'm, in the 
systematic reduction of the E2 transition probabilities for the first excited 2"'" states of Sn 
isotopes, and also in the systematic reduction of E3 transition probabilities for the first 
excited 3~ states of Mo isotopes. It should be noted here that recently the doubly-stretched 
quadrupole pairing interaction was successfully applied also to the microscopic analysis 
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of identical bands in super deformed nuclei, and it was shown that the doubly-stretched 
quadrupole pairing interaction has several advantages compared to the non-stretched and 



stretched ones 41 



In summary, it is clarified that the self-consistent velocity dependent effective interactions 
play crucial roles to recover the local Galilean invariance and eliminate various unphysical 
effects arising from the spurious velocity dependence of the mean potential. For rotating 
nuclei, we can apply similar prescription in order to find the proper effective interactions 
which faithfully take into account the effects of the collective rotation. Results of it will be 
reported in a separate paper. 
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FIGURES 

FIG. 1. El transition strength distribution for a schematic model system of N=Z=20 with 
axially symmetric deformation of /3 = 0.4 calculated in RPA. The continuous strength function, 
representing the strength per unit energy, is given in units of MeV~^ relative to the classical sum 
rule value (CSR). In the model Hamiltonian of Eq.(^2|), the pairing gap is fixed as A = 2.0 MeV 
both for protons and for neutrons. K=0 and K=l modes are shown by the solid and the dashed 
curves, respectively, (a) and (b) correspond to the results obtained with and without the inclusion 
of the dipole-pairing interaction, respectively. 

FIG. 2. El transition strength functions of (a) '^'^^Sm and (b) ^^'^Sm are given in units of 
MeV~^ relative to the classical sum rule value (CSR). The model Hamiltonian is same as that of 
Eq.([85|) except that the dipole-pairing interaction is not included here. In (b), solid and dashed 
curves correspond to the K=0 and K=l modes, respectively. 



FIG. 3. Same as Fig]2|, except that the dipole-pairing interaction is included. 

FIG. 4. Same as Fig.^ for (a) and (b), except that the additional interaction H^^^^ is included. 
For comparison, Lorentzian distributions which fit the experimental data of the photoneutron cross 
section |35] are given in (c) and (d) for ^^^Sm and ^^''Sm, respectively. The experimental Lorentz 
line parameters for the mathematical expression (Ti{E) = _^^y+{ev Y taken from Ref. 



|3| as El = 14.8 ± 0.1, Fi = 5.1 ± 0.2, ai = 339 ± 12 for the best single line fit of ^^^Sm and 
El = 12.35±0.10, Fi = 3.35±0.15, ai = 192±10, E2 = 16.1±0.1, F2 = 5.25±0.20, ^2 = 204±10 
for the best two line fit of ^^^^m, respectively. Here £"s and F's are given in units of MeV while 
it's are given in units of mb. 
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FIG. 5. (a) Excitation energies of the first 2"*" states and (b) E2 transition probabilities in 
Sn isotopes. Results of calculations with and without the inclusion of the quadrupole pairing 
interaction are shown by dot-dashed and dashed lines, respectively. Solid lines correspond to 
the experimental data. In the calculations, pairing gaps are determined from the odd-even mass 
differences both for protons and for neutrons. 

FIG. 6. Same as Fig.|5|, except that the pairing gaps for protons are set to be Ap = in the 
calculations. 

FIG. 7. (a) Excitation energies of the first 3~ states and (b) E3 transition probabilities in Mo 
isotopes. Results of calculations with and without the inclusion of the octupole pairing interaction 
are shown by dot-dashed and dashed lines, respectively. Solid lines correspond to the experimental 
data. 



32 



TABLES 

TABLE L El transition strength distribution for a schematic spherical system with N=Z=20. 
Energies and fractions of the EWSR are shown for some dominant states calculated in RPA. The 
dipole-pairing interaction is included for c and e, but not for a, b, and d. 



= A„ (MeV) 


E (MeV) 


EWSR (%) 


a 0.0 


24.0 


100 


b 1.0 


21.4 
27.5 


8 

41 


c 1.0 


24.0 


100 


d 2.0 


20.8 
23.6 
30.0 
51.3 


16 
24 
78 
13 


e 2.0 


23.9 


99 


TABLE IL Calculated El EWSR values of ^^^'^^^Sm. Strengths integrated over the excitation 
energies from 5.5 MeV to 30.0 MeV are given in units of % relative to the classical sum rule value. 

The model Hamiltonian is given by Eq.(|85|). The columns a and b correspond to the results 

(p) 

obtained without and with the inclusion of the dipole-pairing interaction H^_^, respectively, while 
in c the interaction H^f^ is also included in addition to H^l^. 


Nucleus a 


b 


c 


i^^Sm 116.8 
^^'^Sm (K=0) 41.1 
1545m (K=l) 64.7 


85.9 
31.3 
56.0 


86.2 
30.2 
56.9 
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TABLE III. Energies of the first excited 2+ states and E2 transition probabihties in Sn isotopes 

are given in units of MeV and B{E2)sp, respectively. Results of calculations obtained without and 

(p) 

with the inclusion of the quadrupole-pairing interaction H^J^^ ^"^^ compared with experimental 
data. The pairing gaps adopted in the calculations are taken from the experimental even-odd mass 
differences for a, while those for protons are set to be zero for b. 



a b 





without H^l^ 


with H^l^ 


without i7^^2 


with H^}^ 


Exp. 




N 




B(E2) 


Ei2+) 


B(E2) 


E{2+) 


B(E2) 


E{2+) 


B(E2) 


E(2+) 


B(E2) 


62 


1.72 


18.7 


1.57 


13.3 


1.66 


15.7 


1.57 


10.7 


1.26 


16.2 


64 


1.72 


19.6 


1.57 


14.1 


1.68 


15.9 


1.58 


10.8 


1.30 


15.3 


66 


1.41 


19.8 


1.19 


17.1 


1.37 


15.7 


1.34 


10.6 


1.29 


12.9 


68 


1.23 


23.9 


1.03 


20.4 


1.20 


19.0 


1.21 


12.2 


1.23 


12.7 


70 


1.20 


24.8 


1.00 


21.2 


1.16 


20.0 


1.18 


12.7 


1.17 


11.6 


72 


1.22 


22.7 


1.13 


15.0 


1.17 


18.5 


1.18 


11.9 


1.14 


10.9 


74 


1.35 


18.0 


1.27 


12.9 


1.31 


14.9 


1.27 


10.0 


1.13 


9.3 
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TABLE IV. Energies of the first excited 3 states and E3 transition probabilities in Mo isotopes 

are given in units of MeV and B{ES)sp, respectively. Results of calculations obtained without and 

(p) 

with the inclusion of the octupole-pairing interaction H^J^^ are compared with experimental data. 



The pairini 


5 gaps adopted 


in the calculations are taken 


from the 


experimental 


even-odd mass 


differences. 
















without ij|^3 


with if^ 


3 




Exp. 


N 




B(E3) 




B(E3) 


E{3-) 


B(E3) 


52 


2.74 


51.2 


2.66 


44.3 


2.53 


17 


54 


2.36 


53.4 


2.32 


44.6 


2.23 


24 


56 


1.93 


58.6 


1.97 


47.3 


2.02 


33 


58 


1.91 


60.2 


1.92 


47.8 


1.91 
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